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Abstract: If, as recently reported by the Super-Kamiokande collaboration, the neutrinos are massive, 
the heaviest one uh would not be stable and, though chargeless, could in particular decay into a 



Qh! lighter neutrino ul and a photon by quantum loop effects. The corresponding rate is computed in the 

' standard model with massive Dirac neutrinos as a function of the neutrino masses and mixing angles. 

The lifetime of the decaying neutrino is estimated to be « 10^^ years for a mass « 5 x 10^^ eV. 



If kinematically possible, the vh —>■ i^L e^ e~ mode occurs at tree level and its one-loop radiative 
corrections get enhanced by a large logarithm of the electron mass acting as an infrared cutoff. Thus 
C^ I the vh — ^ t^L e^ e~ decay largely dominates the vh ^ ^l 1 one by several orders of magnitude, 

corresponding to a lifetime ss 10"^ year for a mass ps 1.1 MeV. 
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1 Introduction 

Evidence for the transmutation between the two neutrino species v^ ^-> v^ has been recently reported 
by the Super-Kamiokande collaboration [|l]]. As a consequence, neutrinos could have non-degenerate 
tiny masses, and mixing among different lepton families becomes likely, in analogy with the Cabibbo- 
Kobayashi-Maskawa flavor mixing in the quark sector |^. 

We assume that the neutrino "flavor" eigenstates v^, t^^ and Vr are linear combinations of the three 
neutrino mass eigenstates vi, V2 and z^3 of nonzero and non-degenerate masses mi, m2 and m^ 
respectively according to 



U,\ 



where the 3x3 matrix Uiep is unitary. 

The effective weak interactions of the leptons can now be written as 
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where the charged current Lx is 
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Here I stands for e , /i , r and Vi (with i = 1, 2, 3) are the three neutrino mass eigenstates. 

Although the neutrinos are chargeless, a heavy neutrino vh can decay into a lighter neutrino vl by 
emitting a photon; this decay is entirely due to quantum loop effects. Now, if kinematically possible, 
the mode vr —>■ ^l e"*" e~ largely dominates, because it is governed by a tree diagram and its radiative 
corrections get enhanced, as we will see, by a large logarithm. 

Neutrino oscillation measurements provide constraints usually plotted in the {sm.29ij, 
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m^ I ) plane, where 9ij is one of the three Euler angles of the rotation matrix Ui 



ep- 
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023 ~ 45*^ is suggested by the Super-Kamiokande data and 0i3 ss O'' comes from the CHOOZ data 
[|l|, ||] which give 0i3 < 13°, and also from the Bugey experiment [^, whereas 0i2 is arbitrary. 
Although 012 is likely small ss O", the maximal mixing 0i2 ~ 45*^ may also be possible allowing 
i^e ^^ ^11 (as suggested by the LSND experiment ^ |3|]). 

2 The decay vh -^ ^l1 

The first calculations of radiative neutrino decays have been reported in [||] and [^. 
In the most general renormalizable gauge (conventionally called R^, six Feynman diagrams con- 
tribute to the process lyniP) -^ ^l{p) 7(9). where the photon can be real {q^ = 0) or virtual {q^ / 0); 
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the latter is necessary when we consider the one-loop radiative corrections to u^ 
can be grouped into two sets: four in Figs, la-d: 
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Fig. Ic 



Fig. Id 



and two in Figs. 2a-b: 




AAAAAAAy 




Aaaaaaa; 



Fig. 2a 



Fig. 2b 



Each one is gauge-dependent but it turns out that the ^ dependence cancels out for each group of 

diagrams separately, yielding the overall gauge independence of the physical process. 

We shall give the results in the 't Hooft-Feynman gauge ^ = 1. 

For each diagram, the corresponding amplitude A is written in terms of the effective vertex F^ 

(^] V ULeUh u{p) r^(£) u{P) e*Jq), 
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where the u's are the (Dirac) spinors of the two neutrinos, e* is the photon polarization, e the charge 
of the electron and g the SU{2)l coupling constant. One has Gp/V^ = 5^/8M^. 

The ultraviolet divergences are handled via the procedure of dimensional regularization, going to 
n = 4 — e dimensions. 

The mass m of the lightest (outgoing) neutrino is always neglected, such that the results depend on 
the mass M of the incoming neutrino, the mass M\y of the W gauge boson, and the masses m^ of 
the internal fermions, which will always appear in the dimensionless ratio 
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After expressing the amplitude for each diagram in terms of two-dimensional parametric integrals, 
we restrict ourselves in this section to the case of a real outgoing photon, for which, due to q^e* = 
and to the conservation of the electromagnetic current, only the magnetic form factor proportional to 
if^fiv <f in the effective vertex contributes (see for example [^[§]). The integration over the Feynman 
parameters is made simpler by neglecting M'^ /My^ in the denominators. 



2.1 Computation of the six diagrams 

• Diagram la The corresponding effective vertex r^*^ writes 

where 

Vi (i) =M^[{l-x)+ rix] -M'^xy-q^y{l-x-y) 
and (7 is the Euler constant 7 « .577) 
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In (^ and in the rest of the paper, A is an arbitrary scale coming from the dimensional regularization. 

By translational invariance, T^^ depends only on the four-momentum transfer g'^ and not on P'^; the 
latter may be projected onto the basis formed by three independent four-vectors i a'^'^q^, q^, and 7^^, 
using the following relation valid for ttt, = 0: 
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The X, y integrations for the pure magnetic term yield the contribution of diagram la to the decay 
amplitude vr ^ ''^Ll 



Aia = AoY,UmUl^Fia{^), 
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where 



One gets 
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The singularities of Fia{t) at r^ = 1 are fake: Fia{i) = —5/12 for r^ = 1. 

Formula (|l^) is in agreement with similar calculations [^ for fi^ —>■ e^ j at the limit r^ — > 0, where 
only the linear term in r^ was kept and the logarithmic term neglected. 



If m were not neglected, the M(l + 75) term in (13) would be simply replaced by M(l + 75) + 
m{l — 75). If we keep M'^xy in Vi (i), we will still obtain explicit analytic forms for the F's but the 
results will be complicated and not illuminating. 



Diagram lb Writing in a similar way 
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AAf, = r,<^M(l-y)[(2y-l)P^ + (l-2x-2y)p^] 
The use of ([l0|) transforms the above expression into 
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and, after performing the parametric integration of the purely magnetic term one obtains 
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The singularities of Fii,{£) at r^ = 1 are again only apparent; in fact Fijj{£) = —1/8 for : 
The computations proceed along the same way for the other diagrams. 
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• Diagram Ic 

AAf^ = iMa^'^^ix + y-l)+ <iq^{l - x - y) + {M\x - 1) + ml) 7^ 
gives after the integrations over x and y 
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• Diagram Id 



yields 
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• Diagram 2a Calling 

V2{i) = M^x + mj{l -x)- M^xy - q^ y{l - x - y) 
one has 
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• Diagram 2b 
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2.2 Cancelation of the ultraviolet divergences 

All terms that are ^-independent do not contribute to the amplitude because of the unitarity of Uiep', 
this is in particular the case of the (divergent) terms (2/e + In 47r — 7 — 1/2) in the diagrams la and 
2a. 

The only two remaining divergent diagrams are lb and 2b; however the coefficients of their (£- 
dependent) divergent terms exactly cancel, ensuring the finiteness of the final result. 



2.3 Result for the total amplitude of uh -^ ^lI 



Dropping the constants (-7/12), (-1/4), (-2/3) in Q, @ and @ which, being ^-independent, 
do not contribute to the decay amplitude (see above), we obtain for the sum of the six contributions 

Ylt UHeUl^[Fi.a-dii) + F2.a,bii)] ^e expression 
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where Aq has been defined in (|13[). Our result (g8|) agrees with formula (10.28) for the function /(r) 
in reference [||] (where the three irrelevant constants mentioned above are kept). 
The corresponding decay rate is 
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With the assumptions about Z^^gp and the corresponding mixing angles mentioned in the introduction, 
one finds for M « 5 10"^ eV 
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(30) 



This is to be compared with the experimental lower limit found in [^]. 

The detectability of this decay and its relevance for astronomy has been emphasized for example in 



[101. 



3 The decay vh 
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If kinematically allowed, this decay is governed at tree level by the diagram of Fig. 3, and at the 
one-loop level by ten diagrams: the six previously considered in Figs. 1,2 where the photon, now 
off-mass-shell, decays into an electron-positron pair, and the four box diagrams of Fig. 4 in which the 
W~^ — W~ pair is converted into the e+ — e^ pair. 
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Fig.3 



The tree amplitude 
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can be recast, by a Fierz transformation and using the unitarity oiUiep, into 
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As for the one-loop corrections, a careful examination of all the terms in (|n|), (jT^), (|19|), (^Tj), ( p4| ) 
and (26) for the six vertices r^^_^(£), Fg^ ^(£) shows that the dominant behavior comes from the q^ 
term in (|^) corresponding to Fig. 2a; it exhibits a In r^ -^ oo for r^ ^ contribution, reflecting 
mass singularities (or infrared divergences) of the loop integrals. 

We can track down this divergent behavior by examining the integration limits x = Q and a; = 1 of 
the denominators I?i 2(^). When r^ = 0, an infrared-like divergence occurs if the numerators M2{i) 
lack an x term to cancel the x = integration limit of the xM"^ term in the denominator P2 (^)- This 
happens with the 2y{y - l)q^ term of A/'^^(£) in (||). 

This infrared-like divergence, which arises when there are two massless (r^ = 0) internal fermions in 
the loop, has been noticed a long time ago in the computation of the neutrino charge radius [11]. 
Compared to In r^, all other terms are negligible because they are strongly damped by powers of 
r£, or r" In r^, where n > and r^ < 10"'^. Thus Fig. 2b is damped by r^ In r^, and the four dia- 
grams of Fig. 1 are all strongly damped since an infrared-like divergence cannot occur here: the x = 1 
integration limit of the (1 — x)M^^ in the denominator I?i(^) is systematically canceled by the 
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coming from the integration over the y variable. Explicit x, y integrations of all six ver- 



tices ^ia~.d(^)^^2ab(^) Confirm these features. 

Similar considerations show that the box diagrams of Fig. 4 share the same power suppression 
r" In r£ as the five other diagrams of Figs, la-d and Fig. 2b. The origin of this r^ power suppression 
in all one-loop diagrams except Fig. 2a can be traced back to the fact that they involve two {W, <&) 
propagators; only Fig. 2a and Fig. 2b have one, but the latter nevertheless gets an r^ suppression from 
the $-fermion couplings. 
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Fig. 4 



To summarize, at one-loop, only the 2y{y — l)g^ term in (^) yields an infrared-like divergence 
oc In rg while all other terms get damped by powers of r^. 

The leading g^ Inr^ term of Fig. 2a in the i^h — i^L — 1 vertex cancels the photon propagator 1/q^ in 
Fig. 5 and yields an effective local four-fermion coupling proportional to Gp- The leading contribu- 
tion to the one-loop radiative corrections to the vh -^ vl e^ e~ tree amplitude is accordingly found 
to be 
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which can be put, using again the unitarity of Uiep into a form similar to Atree in (|32|): 
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The sum Atree + Arad = B h uow casy to manipulate when we consider the interference between 
Atree and Arad in |;Sp for the decay rate. 
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From the amplitude B, we compute [ 12] the decay rate Fi = F 
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and find 
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from which one gets 
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where x = m^/M , and G{x), H{x) are the phase-space functions given by 
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To this leading logarithmic radiative correction expressed by ss a In r in ( |36|J38| ), we may also add the 
non-leading (simply a, without In r) electromagnetic connection to the e"*" e~ pair. This non-leading 
correction can be obtained from the one-loop QCD correction to the well known e~^ e~ — > quark-pair 



cross-section, or the t ^ Ur + quark-pair decay rate found in the literature [12]; the only necessary 
change is the substitution Og <-^ 3a/4. Thus, in addition to Fi, we have the non-leading contribution 



the function K{x, x) is tabulated in Table 14.1 of [12]. 



We emphasize that K{x,x) is a spectacular increasing function of x, acting in the opposite direction 
to the decreasing phase-space function G(x). 



The present direct experimental limit on the mass of z/,- is [ ]13| ] m,y^ < 18.2 MeV; if we take, for 
example, the mass of the heavy decaying neutrino to be 1.1 MeV, its lifetime is found to be w 10^^ 
year. 

Other stronger limits (below 1 MeV) mainly come from cosmological arguments [|l4|][15]. 
Finally we note that the virtual weak neutral Z^ boson replacing the virtual photon in Fig. 4 also 
contributes to uh -^ vl e"*" e~. However it can be safely discarded, being strongly damped by 
q^ /M'^ due to the Z^ propagator. 

4 Conclusions 

The recent observation by the Super-Kamiokande collaboration of a clear up-down z/^ asymmetry in 
atmospheric neutrinos is strongly suggestive of v^ — > vx oscillations, where vx may be identified 
with Vr or even possibly a sterile neutrino. These results have many important physical implications. 
In particular, neutrino oscillations mean that neutrinos have a non-vanishing mass, which, according 
to the new data, may be at least as heavy as 5 x 10^^ eV. If a neutrino vh has indeed a mass, it may not 
be stable against decay and could in principle decay into a lighter neutrino, vl, through a cross-family 
electroweak coupling. We have studied two such decay modes, vh ^ vl^ and uh — > v^ e^ e~, and 
found that the latter, which, in contrast to the former, arises at tree level and gets further enhanced by 
large radiative one-loop corrections, is by far the dominant process and may therefore be detectable 
provided that uh has a mass > 2 rrif.. A positive evidence for such decay modes would give a clear 
signal of the onset of 'new physics'. 
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